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Abstract. We present the results of UBVJHKLM photometry of ° spanning the period from 1976 to 2001. 
Studies of the optical light curve have shown no evidence of any stable harmonics in the variations of the stellar 
emission. In the L band we found semi-regular oscillations with the two main periods of ~ 3.3 yr and ~ 11.9 yr 
and t he ful l amplitude of ~ 0™8 and 0?^6, respectively. The colors of the warm dust shell (resolved by Ohnaka 
et al. 2001) are found to be remarkably stable in contrast to its brightness. This indicates that the inner radius 



is a constant, time-independent characteristic of the dust shell. The observed behavior of the IR light curve is 
mainly caused by the variation of the optical thickness of the dust shell within the interval r{V) = 0.2 — 0.4. 
Anticorrelated changes of the optical brightness (in particular with P ft; 3.3 yr) have not been found. Their absence 
suggests that the stellar wind of ° deviates from spherical symmetry. The light curves suggest that the stellar 
wind is variable. The variability of the stellar wind and the creation of dust clouds may be caused by some kind 
of activity on the stellar surface. With some time lag, periods of increased mass-loss cause an increase in the 
dust formation rate at the inner boundary of the extended dust shell and an increase in its IR brightness. We 
have derived the following parameters of the dust shell (at mean brightness) by radiative transfer modeling: inner 
dust shell radius rin ~ llOi?*, temperature Tduat(nn) ~ 860K, dust density Pdust('"in) ft; 1.1 x 10^^" g cm"'', 
optical depth r(V) ft; 0.32 at 0.55 ^m, mean dust formation rate Mdust ft; 3.1 x 10~^MQyr~^, mass-loss rate 
Mgas ft; 2.1 X 10~^ Mq yr~^, size of the amorphous carbon grains ^ 0.01 /im, and B ~ V ~ —0.28. 

Key words, stars; carbon - stars: circumstellar matter - stars: individual: RCrB - infrared: stars 



1. Introduction 

R Coronae Borealis is the prototype of a small group 
of yellow supergiants (about 35 known members in our 
Galaxy) characterized by sudden declines in their optical 
brightness and extremely hydrogen-deficient, carbon- rich 
atmospheres. The visual light curve of ° contains quasi- 
regular low-amplitude variations without any dominat- 
ing periods. These are commonly interpreted in terms of 
stellar pulsations and episodic deep declines, which typ- 
ically last a few months and have amplitudes of up to 
8 magnitudes. Such events are thought to be the result of 
obscuration by dust clouds generated spasmodically and 
blown away from the star by radiation pressure (Clayton 
|1996| and references therein) . 

The near-infrared excess, discovered by Stein et al. 
( |1969D , is assigned to a warm dust shell (blackbody tem- 
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perature of ~ 900 K). This extended dust shell was re- 
solved for the first time by speckle interferometric obser- 
vations with the SAO 6 m telescope (Ohnaka et al. 2001 ). 
The IR excess contributes about 30% of the total flux and 
is permanently present, regardless of the visual brightness 
of the object. The L flux of ° which is mainly due to dust 
emission, varies semi-regular ly with a period of 1260 days 
(Feast et al. 1997 and references therein). The analysis 
of IR AS observations at 60 fim and 100 /im (Gillett et al. 
|l986 ) led to the discovery of a very extended "fossil" shell 
around ; whose diameter and temperature are ~ 18' and 
30 K, respectively. 

In this paper we focus on the investigation of the time- 
dependent characteristics of the extended warm dust shell 
and their interpretation. In the next section we present 
the results of our photometric observations of ° ([/ BV in 
1994-1999 and JHKLM in 1983-2001) and the complete 
tables of photometric observations of the star in 1976- 
2001. The basic characteristics of the optical and IR ra- 
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diation of the star during its bright state are analyzed 
in Sect. 3. Radiative transfer modeUng of the warm dust 
shell at its average brightness is described in Sect. 4. The 
variations of the dust shell brightness and colors, and, cor- 
respondingly, its structural parameters, are investigated in 
Sect. 5. Our results are discussed in Sect. 6. 

2. Observations 

Our U BV observations of ° were carried out in the pe- 
riod 1994-1999 with the 0.25 m telescope of the Fairborn 
observatory (Genet et al. 



Table 3. The average magnitudes of ° in bright state. 
< m >, < m >L,max, and < m >L,min denote the mag- 
nitudes during the average, maximum and minimum star 
brightness in the L band 



1987). The results of these 



observations together with those of earlier observations 
(1985-1993) by Fernie & Seager (1994) are summarized 
in Table 1 in http://infra.sai.msu.ru/ftp/rcrb. This table 
contains 1170 measurements with the accuracy of 0'."02 
or better. The standard star was HD 141352: V = 7.476, 
B-V = 0.439 and C/ - B = -0.003. 

JHKLM observations were carried out between 1983 
and 2001 with the 1.25 m telescope of the Crimean 
Station of the Sternberg Astro nomic al Institute (for tech- 
nical details see Nadjip et al. 1986 ). The results our ob- 
servations, togeth er wi th those obtained in 1976-1979 
(Shenavrin et al. 1979 ), are sum marized in Table 2 in 
http : //infra.sai. msu.ru/ftp/rcrb . This table contains 252 
measurements with the accuracy of 0™03 or better. The 
standard star was BS 5947: J = 2.09, H = 1.60, K = 1.30, 
L = 1.12 and M = 1.35. 

The VJKLM light curves of ° during the period 
1983-2001 are presented in Fig.|l|. Also plotted are a 
few estimates in the V band taken from the internet 
page of the association of amateur astronomers VSOLJ 
(http://www.kusastro.kyoto-u.ac.jp/vsnet) and the mea- 
surements reported by Feast et al. ( 1997| ) in the J and 
L bands. As one can see in this figure, a number of deep 
declines of the star's visual brightness were observed in 
1983-1989 and 1995-2001, while during the interval 1989- 
1995 the V flux was almost constant. The very different 
behavior of the L and M light curves will be discussed in 
subsequent sections. 

3. Basic characteristics of R CrB in the bright 
state 

We call the state of ° as bright, if < 6.05 and J < 5.25. 
In this state the contribution of the dust to the optical ra- 
diation of ° is negligible. The average UBV JHKLM and 
RI magnitudes of ° in the bright state are summarized 
m Table|. These values were derived from Tables 1 and 2, 
and th e data presented by Fernie ( 1982| ) and Fernie et al. 
(119861) . The maximum magnitudes of ° during the whole 
observing period are V « 5.69 and J « 4.89. 

3.1. Periodicity 

The Fourier analysis of the visual hght curve (1985-1999) 
during the bright state has revealed a number of harmon- 
ics with the full amplitude of ~ 0™05 concentrated in the 
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interval of 30-60 days. The maximum amplitude of 0^"07 
has the mode with the period P « 39.8 days (see Fig.|). 
The corresponding harmonic with the full amplitude of 
~ 0™025 was also found from the analysis of the B — V 
curve. The brightness maxima follow the color minima 
with a lag of 3 days, and the ratio of the amplitudes of 
{B — V) and V harmonics is about 0.3. No trend has been 
found in the optical brightness. It should be noted that 
on a time-scale of several hundreds of days ° can occa- 
sionally exhibit well-developed pulsations with periods of 
35, 44 and 51 days (Fernie & Seager 1994 and references 
therein) . 



The low-amplitude quasi-regular variability in the J 
band resembles that in the optical region. A period search 
has revealed three main harmonics with periods of 43.7, 
33.2, and 28.1 days and an am plitude of about 0™08. At 
the same time Rao et al. ( 1999 ) have found that the radial 
velocity curve exhibits the only one stable pulsation period 
of about 42.7 days. 

3.2. Spectral energy distribution 

The spectral energy distribution (SED) of ° relating to its 
average brightness in bright state (see Table ^ is shown 
in Fig.|. 

The stellar magnitudes were corrected for the interstel- 
lar reddening accord ing t o the color excess E{B — V) — 
0.05 (Asplund et al. |l997| ). The flux at 0.23 and 0.25 /^m 
was measured with lUE on 31.05.1991 (see INES data 
from the lUE satellite), when the star brightness in the 
visual and near IR was close to its average values in bright 
state presented in Table||. 

To reconstruct the SED beyond 10 /zm we used obser- 
vation of ° with IRAS on 12-13.09.1983 (Rao & Nandy 
1986). The L flux of the star on these dates was smaller 
than its average value {AL ~ 0.25). To account for this 
difference the fluxes measured by IRAS at 12, 25, 60 and 
100 /im were increased by 8%, taking into account that 
A/(3.5Aim)/A/(ll^m) « 3 (Forrest et al. p^72|). 
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Fig. 1. VJKLM light curves of ; The dots represent the observations reported in this paper. The crosses in the visual 
light curve denote the data of VSO LJ (http://www.kusastro.kyoto-u.ac.jp/vsnet), and in the J and L light curves, the 
data reported by Feast et al. (1997). The vertical bars in the L and M panels mark the times of the minima calculated 
from Eqs. ( p4[ ) and ([T5|), respectively. 



Integrating the SED presented in Fig.^ we evaluate 
the average bolometric flux of ° during its bright state 
as -Fboi.mean ~ 1-55 X 10~ erg cm~^ s" ^ About 85% of 
the radiation detected from the star is emitted within the 



interval 0.36 /im< A ^ 5/im, which is completely cov- 
ered by our photometric observations. The contributions 
at A ^ 0.36 Aim, A > S/xm and A <; 12 ^m are 4%, 10% 
and 3%, respectively. In particular, this indicates that the 
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Fig. 2. The V and L Hght curves of ° folded with the 
periods 39.82, 1206 and 4342 days, from bottom to top 



possible errors in the correction of the IRAS data cannot 
significantly change the estimates of the bolometric flux. 

3.3. Infrared excess 

The strong infrared excess (i.e. the bolometric flux of the 
dust shell) can be evaluated as i^dust — ^boi — -F'*,boi- Here 
i^boi is the observed bolometric flux and F*,boi is the bolo- 
metric flux emitted by the central star and passed through 
the dust shell. 

To evaluate the infrared excess we have to reconstruct 
the observed SED of the star (F*(A)). The upper limit of 
the contribution of the dust shell in the J band can be de- 
rived directly from our observations in 1996. In this year a 
deep decline was observed (A J > 3.6), while its L bright- 
ness was close to the average value presented in Table ||. 
Attributing the observed J flux in minimum to the dust 
shell we conclude that its contribution to the total flux out 
of minimum is smaller than 4%. This indicates that within 
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Fig. 3. SED of ° for average brightness in bright state 
(circles). Long- and short-dashed lines denote the SED of 
the star and the dust shell, respectively. The solid line is 
the fitted SED (see text) 



the interval A < 1.25 /im the total flux is approximately 
equal to F*(A). 

To evaluate F*(A) in HKLM we use the intrinsic col- 
ors of ° stars calculated by Asplund et al. ( 1997 ). They 
depend on the effective temperature and the gravitational 
acceleration on the star surface. For ° we take Toff ~ 



6750 K, L 



bol 



10*^0, \ogg w 0.5, radius i?, « 73.4i?, 



■0' 



and correspondingly the intrinsic colors {H — K)q « 0.07, 
{K - L)o ~ 0. 06 an d (L - M)o ~ 0.03 (Asplund et al. 
1997| ; Asplund ^OOOl ). These colors were corrected for the 



interstellar and circumstellar reddening although both of 
these corrections have only a minor effect on the evalua- 
tion of F*^boi- For A > 4.8 /Ltm we approximate F^{X) by 
the blackbody spectrum with the effective temperature 
6750 K. 

The derived SED of the star is presented in Fig.|^ 
(long-dashed line). Integrating F^{X) we find i^*.boi ~ 

1.17 X 10"'^CrgCm-2g-l ^j^^J Fdust,mcan ~ 3.8 X 

10^^ ergcm^^ s^^. This means that about 24% of the stel- 
lar radiation is absorbed and re-emitted by the dust shell. 
Assuming the geometry of the dust shell to be spheri- 
cally symmetric, one can estimate its optical depth as 
Tcff.mcan ~ 0.28. For ycUow stars and amorphous carbon 
grains the value of the Tcff.mcan differs from the optical 
thickness at 0.55 /im (t(V)) by a few percent. 
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V - J = 0.418F~ 1.62. 



(5) 



Fig. 4. {U~B,V), {B-V,V), {V-R,V), (V-I^V) and 
{V — J,V) diagrams. The soHd lines denote the approxi- 
mations by linear equations (|l|-||) 



3.4. Color-brightness diagrams of the star 

The derived color-brightness diagrams (U ~ B,V), {B ~ 
y, V),{V- R, V),{V- 1, V) and {V ~ J, V) are presented 
in Fig. ^. These diagrams have been constructed using the 
data from Tables 1 and 2 and the values reported by Fernie 
( |1982D and Fernie et al. ( [19861 ). For the diagram {V~J, V) 
we have chosen the estimates which were obtained in all 
three bands on the same date. Only the data relating to 
the star brightness in the L band smaller than the average 
value were used. The J flux has been reduced by 4% in 
order to exclude the contribution of the dust shell. The 
linear approximation of the diagrams yields 



[/ - B 0.402t/ - 2.31, 
B-V ^ 0.280X^-1.08, 
F-i? = 0.301^-1.36, 
1/-/ = 0.430^-1.94, 



(1) 
(2) 
(3) 
(4) 



where 5.79 < V < 6.05. It follows from these equations 
that the star reddens as its light declines. 



4. Modeling of the warm dust shell at its mean 
brightness in bright state 



Using the quasi-diffusion method of Leung ( 1976), as im- 
plemented in the CSDUST3 code (Egan et al. |l988| ), we 
have constructed radiative transfer models which match 
the SED of e The input parameters of the models were the 
same as described by Yudin et al. ( 2001 ) . 

As a first approximation we have assumed that the 
warm dust shell is spherical and consists of amorphous 
carbon grains with the radius a^^ = 0.01 /im. The effi- 
ciencies of absorption and scattering of the dust grains 
were calculated using the Mie theory. The refractive in- 
dex was taken from Jager et al. ( 1998| , sample cellOOO). 
The intrinsic SED of the star emission, F^:fi{X), was cal- 
culated as F,_o(A) = F^{X)e'^'^^\ where t(A) is the total 
optical depth at a given reference wavelength. The refer- 
ence wavelength was chosen to be 0.55 /im, and the model 
SED is normalized to the observed one in the V band. 

We have derived the following parameters of the dust 
shell which refer to its mean brightness in bright state. 
The optical depth at 0.55 /im and 2.2 /im are t{V) w 0.32 
and t{K) w 0.05, respectively. The inner radius of the 
dust shell is nn ~ 110 i?*. The temperature and the 
dust density at the inner radius of the dust shell are 
rdust(ri„) « 860 K and /5dust(nn) » 1-1 x lO^^Og^j^-a^ 
respectively, and the density distribution of the dust shell 
is p{r) oc r~^. 

Using these parameters we have evaluated the dust 
formation rate at the inner boundary of the dust shell as 
-^j'^dust ~ 3.1 X 10~^ Mq yr~^. If one assumes that almost 
all carbon is bound in the dust grains and uses the chemi- 
cal com position of the atmosphere of ° from Asplund et al. 
( ^OOOl) , the gas density at the inner radius of the dust shell 
is Pgas(nn) ~ 7.3 X lO^^^gcm^^. Correspondingly, with a 
stellar wind velocity of Vcnv ~ 45kms~^, which is typical 
for stars of this type, one finds a mass outflow rate in 
of Mgas ~ 2.1 X lO^^MQyr^^. If, however, a fraction of 
the carbon and almost all of the oxygen are bound in the 
CO molecule, the corresponding estimates of /?gas(nn) and 
Afgas turn out to be larger by a factor of 2.5. 

As follows from our modeling, the visual brightness 
of the dust shell is F « 14"^6 and its color indices are 
U - B w -0.65, B -V K -0.28, F - i? « -0.41 and 
i? - / « 1.47. The SED of the dust shell is plotted in 
Fig. H (short-dashed line). The radiation of the dust shell 
in U BVR is dominated by scattering. Since the radius of 
the grains is significantly smaller than the wavelength, the 
scattering is governed by the Rayleigh law (csca oc A~^). 
This explains the blue colors of the dust shell in this part 
of the spectrum. In the / band the thermal emission of 
the grains begins to dominate over the scattering. 
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The visual magnitude of ° observed during the deepest 
dechnes is about 14™. Thus, the visual magnitude of the 
dust shell should be fainter than 14™. This was the reason 
for choosing a small radius of the dust grains 0.01 /^m) 
in our calculations. If one assumes that the radius of the 
dust grains is larger, the visual brightness of the dust shell 
becomes greater than that of the star brightness in deep 
minima. For instance, increasing the radius of the dust 
grains by a factor of 2.5 (i.e. to 0.025 /xm), one finds a 
visual magnitude of the dust shell of F « 11™5. This 
reflects the fact that the albedo of dust grains rapidly 
increases with their radius. 

We have also found the colors of the dust shell to be 
sensitive to the value of rin- Variations of this parameter of 
only 10% lead to variations of the IR excess color indices 
{H - iir)dust and [K - L)dust by ~ IJ^I. At the same 
time these indices are observed to be remarkably stable in 
contrast to the brightness of the dust shell (see Sect. 5). 
In the context of our model this means that the value of 
the inner radius rin is a fairly constant, time- independent 
characteristic of the extended warm dust shell. It should 
be noted, however, that the absolute value of rin depends 
on the optical properties of the dust grains included in the 
model calculations. If we used the extinction of amorphous 
carbon obtained by Colangeli et al. ( 199^ ) (AGAR sample) 
instead of cellOOO (Jager et al. 199S), the value of rin would 
be smaller by a factor of ~ 1.4. 

Finally, according to Reimers equation (Reimers 1975|) , 
the average mass loss by a yellow supergiant can be esti- 
mated as 



M = 4 X 10 



-13 



M0yr-\ 



(6) 



where L*, _R* and M, are the luminosity, the radius and 
the mass of the star expressed in the solar units, and the 
dimensionless parameter is limited to 0.3 < < 3. 
Using the corresponding parameters of ° we find M « 
4.2 X 10~^?7Mq yr~^. This is in agreement with the mass 
outflow rate derived from our model provided rj « 0.5. 
Therefore, the mass loss rate of ° agrees with the Reimers 
equation and no superwind is required for the interpreta- 
tion of the dust shell formation. 

The large value of nn ~ 110 i?,, which does not change 
significantly with time, suggests that the extended warm 
dust shell of ° originates in a similar way as in the case of 
ordinary supergiants, i.e. it is caused by the stellar wind 
and at a large distance from the star. In addition to this 
common type of dust shell of supergiants, the circunistel- 
lar shell of ° contains dust clouds which condense some- 
where in the vicinity of the star from high density gas 
clumps (Feast 1997). These additional dust clouds can ex- 
plain the deep minima in the visual light curve. If the 
dust clouds are considered as an additional source of the 
IR emission, then the effective size of this source is much 
smaller then the size of the above extended dust shell, 
which condenses from the stellar wind. Until now all ef- 
forts to pick out the emission of this consortium of the dust 



Table 4. The bolometric flux, the infrared excess and the 
optical thickness of the dust shell of ° during maximum 
and minimum of L brightness 



erg cm" 



-fdust 



Toff 



erg cm 



Maximum 
Minimum 



1.90 X lO"'' 
1.24 X lO"'' 



6.1 X 10-** 0.39 
2.1 X 10"** 0.19 



clouds from the total IR excess were unsuccessful (Feast 



1997) 



5. Variability of the warm dust shell 

5.1. The dust shell during maximum and minimum 
brightness 

The UBVRIJHKLM magnitudes of ° during its max- 
imum (in 1994 April-May) and minimum (in 1999 May- 
June) brightness in L are presented in the second and 
third columns of Table ^, respectively. The corresponding 
bolometric flux and infrared excess of ; and the optical 
thickness of the dust shell derived using these data are 
presented in Table ^. 

As follows from Table ^ the amplitude of variations of 
the optical thickness of the dust shell in ° is Ar « 0.2. At 
the same time, the value of the color index (K — L)dust of 
the IR excess, estimated during maximum and minimum 
brightness in the L band, proves to be the same: (A' — 
i)dust ~ 2"^35. This indicates that the variations of the 
optical depth of the dust shell cannot be explained by the 
variations of its inner radius, but reflect the changes of the 
dust formation rate at a certain distance from the star. 



5.2. Colors and the optical thickness of the dust shell 

Fig. H shows the color diagrams (i?dust, ^dust), 
(K dust, L dust), and (Afdust, ^dust)- The dots in these 
diagrams denote IR observations of ° when the star was 
in its bright state. It was assumed that the contribution 
of the dust shell to the flux in J band is negligibly small. 
The stellar magnitudes in HKLM and the infrared 
excess have been derived using the method described 
above. The linear approximation of the diagrams gives 



-f^dust — 1.22L 



dust 



4.22, 



Kdust = 0.97Ldust + 2.35, 
Afdust = 0.98Ldust - 0.66. 



(7) 
(8) 
(9) 



For the mean L brightness, -/jdust.mcan ~ 2.49, we find {H— 

^)dust,mcan = 4.76 ± 0.4, {K - L)dust,mcan = 2.28 ± 0.15, 

and (L - M)dust .mean — 

0.70 ±0.12. 

As follows from Eqs. (0-||), the variations of the dust 
shell brightness are not accompanied by remarkable vari- 
ations of its colors. Hence, there is a proportionality be- 
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Fig. 5. The (ffdust, ^dust), (^''dust, ^dust) and 
(Afdust, ^dust) diagrams of the emission of the dust 
shell. The solid lines denote the approximations by linear 
equations 



tween the bolometric stellar magnitude of the dust shell, 
JTidust.boi, and its L magnitude: 



'Tldust.bol — 



dust 



4'.^55. 



(10) 



We can calculate the function t^s = /(i^dust) from 
the expression i^dust = ^boi(l — ea^p(— Toff)). Finally, us- 
ing a quadratic polynomial approximation of the functions 
idust = f{L) for the data of our IR observations, we find: 



mdust.boi ^ 0.046L^ + 0.88L + 4.68, 
Toff ~ 0.176^2 - 1.20L + 2.14. 



(11) 
(12) 



The function Toff (<) for the period 1994-2001 is plotted in 
Fig.^. Our modeling has shown that such a sharp transi- 
tion from the decrease to the increase of the optical depth 
requires an abrupt change of the dust formation rate close 
to the minimum. This means that there was a sharp den- 
sity enhancement in the stellar wind. The solid line in 



Fig. H shows the function roft(i), calculated with the as- 
sumption that at minimum the rate of dust formation 
sharply increases in the minimum of Teff by a factor of 
4.3, and then rises oc t^-^ during about 560 days until the 
maximum value of the optical depth is reached. 

The dust formation rate depends on the gas density 
at the inner boundary of the dust shell and, therefore, on 
the intensity of the mass-loss. The time lag between the 
increase of the stellar wind and the increase of gas density 
at the inner radius of the shell is 



VllOR* 



45 kms 



yr. 



(13) 



Interestingly, this 4- year wind travel time can probably 
be seen in the V and L light curves. The comparison of 
the V and L light curves shows that the period 1994-1999 
of decreasing L and M brightness follows approximately 
4 years after the period of constant V magnitude between 
1989 and 1995. We suggest the following explanation for 
this. We assume that during phases of increased stellar 
activity there is a simultaneous increase of both (l)the 
stellar (supergiant) wind and (2) the rate of the forma- 
tion of (eclipsing) dust clouds (and we assume that the 
eclipsing dust clouds cause the deep minima in the V light 
curve). If these two effects are correlated, then 1983-1989 
(with several V minima) was a period of increased activ- 
ity and, therefore, increased stellar wind. After the wind 
travel time of approximately 4 years, this increased stellar 
wind arrived at the extended dust shell region and then 
caused the L band brightness increase. On the other hand, 
in the period 1993-1999 the L band brightness was de- 
creasing because in this period only lower intensity stellar 
wind arrived, which was produced during the low-activity 
period 1989-1995 (high V brightness). 

5.3. Periodicity 

In order to search for the periodicity in the variations of 
the dust shell brightness, i.e. of the mass-loss rate, the L 
light curve has been analyzed using the code developed by 
Yu.K. Kolpakov (http://infra.sai.msu.ru/prog/kolpakov). 
We combined our observations with the data presented 



by Str ecker ( [1975| ; 35 points during 1968-1974) and Feast 
et al. ( |1997| 24 points during 1983-1991), and removed 
the points obtained during visual declines. Two harmon- 
ics of semi-regular variations of the star brightness in the 
L band were derived: Pi ~ 1206 days, ALi « 0^'^8 and 
P2 ~ 4342 days, AL2 « 0^6. The corresponding phase 
diagrams are presented in Fig.|^. The dates of minima of 
these harmonics are given by the following expressions: 



min(Li) = 2442986 + 1206'^E, 
min(L2) = 2442860 + 4342'^E. 



(14) 
(15) 



These dates are marked by vertical bars on Fig.|l|. The 
minima of these harmonics coincided in 1999, when the 
star brightness in L and M was in unusually deep min- 
imum. It should be noted that the existence of the two 



8 



B. Yudin et al.: U BV J HKLM photometry and modeling of R Coronae Borealis 



0.4 



0.3 



0.2 





1 


i 1 




«• 
t 
















• • 






9** 


t 

• 














» 






* 










• 




••• 


\ 3 








f \ 9 








V / 








\ >• 








/• 








^ / 








M / 

V 











1 



10000 



11000 
JD (+2440000) 



12000 



Fig. 6. Optical thickness Toff of the dust shell versus time 
(see text) 



dominant modes superimposed with irregular variations 
makes the L light curve quite complex. Moreover, the 
shape of the harmonics differs from the simple sine wave 
shown in Fig. ^ For example, the rise of the L brightness 
in 1984 and 1988 was rather sharp compared to the subse- 
quent declines (Fig.|l|). The values of the derived periods 
are very close to those of the mean time interval between 
the i ndividual visual declines of ° ~ 1100 days (Jurcsik 
199^) , and between groups of visual minima, ~ 4400 days 



(Rosenbush ^00l[) 



In the bright state of ° the L variability is mainly 
caused by changes of the optical depth of the shell. Hence, 
if the dust shell is spherical, we would expect to find the 
corresponding semi-regular variations in the visual star 
brightness with period Pi « 1206 days, which are anti- 
phased to those in the L band. According to Eq. ( p^ 
their full amplitude would be AV « 0.3 mag. However, 
a period search has not revealed this harmonic in the V 
light curve. This indicates that the density of the stellar 
wind (and, correspondingly, the shell density) along the 
line of sight is smaller than the average one, i.e. the dust 
shell is asymmetric. 

In this case the derived optical depth of the dust shell 
can be considered as the sphere-averaged one. If the dust 
shell is optically thin, then its bolometric flux is propor- 
tional to the optical depth averaged over the sphere. Since 
the shell's SED is primarily defined by the value of the 
inner radius (i.e., T'dust(?'in)) our above results concern- 
ing the inner radius (its size and time-independence) and 
mass-loss rate would be valid. 



6. Concluding remarks 

We have presented optical and IR long-term monitoring 
and radiative transfer modeling of : These studies have al- 
lowed the derivation of various time-dependent properties 
of both (i) the star itself at bright state and (ii) an ex- 
tended dust shell (similar to that observed around many 
supergiants) which condenses from the supergiant wind 
at a large distance of approximately 100 stellar radii from 
the star. This extended dust shell is larger than the sug- 
gested region of dust clouds (distance approx. 2-30 stellar 
radii, see, e.g. Clayton [199^ , Fadeyev [l98^ ) , which causes 
the deep minima in the visual light curve. The extended 
dust shell (radius ~ 19 mas) was first resolved by speckle 
interferometric observations with the SAO 6 m telescope 



(Ohnaka et al. |200lD . 

In the V and L light curves probably a ~ 4- year wind 
travel time from the stellar surface to the extended dust 
shell with a radius of approximately 110 stellar radii can 
be seen. The comparison of the V and L light curves shows 
that the period 1994-1998 of an increased (and decreas- 
ing) L and M brightness started approximately 4 years 
after the end (in 1990) of a period with many deep min- 
ima in the visual light curve. This can be explained if we 
assume that during phases of increased stella r magnetic 
activity (see studies by Soker & Clayton 1999) there is a 
simultaneous increase of both (i) the supergiant wind and 
(ii) the rate of formation of dust clouds. If this assumption 
is true, the increased stellar wind period ended in 1990. 
The L band brightness increased until 1994. This can be 
explained by the increased stellar wind until 1990 and a 
wind travel time of 4 years from the stellar surface to the 
extended dust shell region. In the period 1994-1999 the L 
band brightness was decreasing because in this period only 
lower intensity stellar wind arrived, which was produced 
during the low-activity period 1990-1995. We believe that 
the variation of the stellar wind can be considered as an 
indirect argument in favor of active phenomena on the 
surface of ; 
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